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We examine how dark energy constraints from current observational data depend on the analysis 
methods used: the analysis of Type la supernovae (SNe la), and that of galaxy clustering data. 
We generalize the flux-averaging analysis method of SNe la to allow correlated errors of SNe la, in 
order to reduce the systematic bias due to weak lensing of SNe la. We find that flux-averaging leads 
to larger errors on dark energy and cosmological parameters if only SN la data are used. When 
SN la data (the latest compilation by the SNLS team) are combined with WMAP 7 year results 
(in terms of our Gaussian fits to the probability distributions of the CMB shift parameters), the 
latest Hubble constant [Hq) measurement using the Hubble Space Telescope (EST), and gamma ray 
burst (GRB) data, fiux-averaging of SNe la increases the concordance with other data, and leads 
to significantly tighter constraints on the dark energy density &t z = 1, and the cosmic curvature 
f2fe. The galaxy clustering measurements of H{z — 0.35)rs{zd) and rs{zd)/ Da{z = 0.35) (where 
H{z) is the Hubble parameter, Da{z) is the angular diameter distance, and rs{zd) is the sound 
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combined. Current data are fully consistent with a cosmological constant and a fiat universe. 
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I. INTRODUCTION 

Solving the mystery of cosmic acceleration [l], [2[ is one 
of the most important challenges in cosmology today. 
Current observational data are not sufficient for differ- 
entiating two likely explanations for the observed cosmic 
acceleration: dark energy, and the modification of gen- 
eral relativity. For recent reviews, see 

There are a number of powerful direct probes of dark 
energy: Type la supernovae (SNe la) P, Q; galaxy clus- 
tering (GC). especially the baryon acoustic oscillations 
(BAG) [H [11; and weak lensing of galaxies [H 
These methods are complementary to each other, and 
each method has its own systematic uncertainties. Gther 
data, e.g., gamma ray bursts can help strengthen 

the dark energy constraints. 

The cosmic microwave background (CMB) anisotropy 
data provide strong priors on cosmological parameters 
(see, e.g., [ll|)- Direct measurements on the Hubble 
constant (see, e.g., [l^) also help break the degeneracy 
amongst the dark energy and cosmological parameters. 

Much progress has been made since the discovery of 
cosmic acceleration in 1998. However, current data are 
still rather limited in constraining dark energy. In par- 
ticular, different analysis methods of the same data could 
lead to significantly different constraints on dark energy. 
This is probably due to the fact that different methods 
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have different residual systematic biases, and they of- 
ten assume different priors on cosmological parameters 
as well. 

In this paper we examine these issues by studying how 
dark energy constraints from current observational data 
depend on the analysis methods used. In particular, we 
explore the impact on the overall dark energy constraints 
of (1) using flux-averaging to minimize the weak lensing 
bias in the analysis of SNe la; (2) using radial and trans- 
verse distance constraints derived from GC data without 
assuming CMB priors. 

We describe our method in Sec. II, present our results 
in Sec. HI, and conclude in Sec. IV. 



II. METHOD 

We only use methods that give geometric constraints 
on dark energy. The constraints on the growth rate of 
cosmic large scale structure are degenerate with the ge- 
ometric constraints (see, e.g., [13, IIH); current data do 
not allow the determination of both without strong as- 
sumptions, e.g., assuming that general relativity is not 
modified. 

Geometric constraints on dark energy are derived from 
the measurement of distances. The comoving distance to 
an object at redshift z is given by: 

r{z) = cH^' \nk\-'/hmn[\n,\'/^r{z)], (1) 

n-A = [4^y Eiz)^Hiz)/Ho 
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where sinn(a;) = sin(x), x, sinh(x) for fl^ < 0, fifc = 0, 
and f2fc > respectively; and the expansion rate of the 
universe H{z) (i.e., the Hubble parameter) is given by 



where z and ^hci are the CMB restframe and heliocentric 
redshifts of the SN. 

For a set of N SNe with correlated errors, we have [23 



/aV = Am^ • • Am (7) 

H'{z)^(^] (2) 

^ ^ where Am is a vector with N components, and C is the 

+ z)^ + nr{l + zf + r2fc(l + zf + ^xX{z)\ ,N X N covariance matrix of the SNe la. 

Note that Am is equivalent to A/iq, since 



where Qm + fir + + = 1, and the dark energy 
density function X{z) is defined as 



X{z) 



Pxjz) 
PxiO)- 



(3) 



Note that fi^ ^ , thus the il^ term is usually omitted 
in dark energy studies, since dark energy should only be 
important at late times. 



A. Analysis of SN la Data 

SN la data give measurements of the luminosity dis- 
tance dL{z) through that of the distance modulus of each 
SN: 



fJ-o 



Af = 5 log 



dLjz) 
Mpc 



25, 



(4) 



where m and M represent the apparent and absolute 
magnitude of a SN. The luminosity distance (Il^z) = 
(1 + z)r(z), with the comoving distance r(z) given by 
Eq.©. 

We use the compilation of SN la data by Conley et 
al. (2011) 22], which include the SNe la from the first 
three years of the Supernova Legacy Survey (SNLS3), 
the largest homogeneous SN la data set. We compare 
two methods for using these SNe in constraining dark 
energy: with flux-averaging of SNe la, and without. 



1. SN la Data 

For a set of 472 SNe la, Conley et al. (2011) [13 
give the apparent B magnitude, tob, and the covariance 



matrix for Am = tub 



with 



mmod = bXog^Q Vl{z\s) - a{s - \) + pC + M, (5) 

where I?i(z|s) is the luminosity distance multiplied by 
Ho for a given set of cosmological parameters {s}, s is 
the stretch measure of the SN light curve shape, and C is 
the color measure for the SN. is a nuisance parameter 
representing some combination of the absolute magni- 
tude of a fiducial SN la, M , and the Hubble constant 
Hq. Since the time dilation part of the observed lumi- 
nosity distance depends on the total redshift Zhoi (special 
relativistic plus cosmological), we have [1^ 



'Dl{z\s)=c '^Hq{1 + ZY,ci)r{z\s), 



(6) 



A?7i — TUB — TOmod = [inB + a{s — 1) — /3C] — M. (8) 



The total covariance matrix is [22 



(9) 



with the diagonal part of the statistical uncertainty given 

by m 



D 



stat,ii 



^ int. 



Hi + z,) 



_z,{l + Zi/2)lnl0_ 



Icnsing ' *^host correction 
2 

<7Z , + a~aZ 



2apCsc,^, (10) 



where CmBs.i, CmgCi, and Csc.i are the covariances be- 
tween niB, s, and C for the i-th SN. Note also that 
(rf j includes a peculiar velocity residual of 0.0005 (i.e., 
15bkm/s) added in quadrature [2^ . 

The statistical and systematic covariance matrices, 
Cstat and Csys, are generally not diagonal [13], and are 
given in the form: 



-'sys 



= Vr 



-l3^Vb + 2aVoa- 2pVob - 2apVab ■ 

(11) 

where Vq, K, Vb, Vba, Vbb, and Kb are ma- 
trices given by the SNLS data archive at 
https:/ /tspace. library, utoronto.ca/handle/ 1807/24512/. 
Cstat includes the uncertainty in the SN model. Cgys 
includes the uncertainty in the zero point. Note that 
Cstat and Csys do not depend on M, since the relative 
distance moduli are independent of the value of M. [13 • 
We refer the reader to Conley et al. (2011) [13 for 
a detailed discussion of the origins of the statistical and 
systematic errors. As an example, we note that the cor- 
relation of errors on different SNe arises from a statistical 
uncertainty in the zero point of one passband, e.g., tm- 
This directly affects all SNe with tm measurements due 
to K-corrections (restframe B to tm), and indirectly af- 
fects even the SNe without tm measurements through 
the empirical SN models by changing the templates and 
the measured color-luminosity relationship. 



2. The Recipe for Flux- Averaging of SNe la 

Because we live in a lumpy universe, SN la observations 
can be misinterpreted if the effect of gravitational lens- 
ing is not properly accounted for (see, e.g., [13|)- Flux- 
averaging of SNe la was proposed to reduce the effect 



3 



of the weak lensing of SNe la on cosmological parameter where 



estimation [25|. The basic idea is that because of flux 
conservation in gravitational lensing, the average magni- 
flcation of a large number of SNe la at the same redshift 
should be unity. Thus averaging the observed flux from a 
large number of SNe la at the same redshift can recover 
the unlensed brightness of the SNe la at that redshift. 

Wang & Mukherjee (2004) ^ and Wang (2005) g3 
developed a consistent framework for flux-averaging SNe 
la. Wang & Mukherjee (2004) [2^ gave the recipe for 
flux-averaging SNe la in the absence of correlated errors. 
Here we modify and generalize the recipe to allow corre- 
lated errors of SNe la. 

As described in [2^ , the fluxes of SNe la in a redshift 
bin should only be averaged after removing their redshift 
dependence, which is a model-dependent process. For 
statistics using MCMC or a grid of parameters, here are 
the steps in flux-averaging: 

(1) Convert the distance modulus of SNe la into 
"fluxes", 
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-(/xr''(^<)-25)/2.5 ^ 



^data 



Mpc 



(12) 



(2) For a given set of cosmological parameters {s}, ob- 
tain "absolute luminosities", {£(z/)}, by removing the 
redshift dependence of the "fluxes", i.e.. 



C{zi) = dl{zi\s)F{zi). 



(13) 



(3) Flux-average the "absolute luminosities" {C}} in 
each redshift bin i to obtain I'C' j: 



' 1=1 



1 



— Vz^^) 

' 1=1 



'■I 



(14) 



(4) Place C at the mean redshift Zi of the i-th redshift 
bin, now the binned flux is 



F{z,)=C'/dl(z,\s) 



(15) 



(5) Compute the covariance matrix oi F(zi) and F{z 



Cov[F(z,),F( 



(16) 



In 10/2.5 



dLizi\s)dLizj\s) 



^ ^ /:(z«)£(zH))(AMrn4'0AMg^'^(^^^)) 

1=1 m=l 

where (A^Q'''''(z}''')A/io''*'*(zm^)) is the covariance of the 
measured distance moduli of the l-th SN la in the i-th 
redshift bin, and the m-th SN la in the j-th redshift bin. 
C{z) is deflned by Eqs.JISl) and (fT3| . 

(6) For the flux-averaged data, |F(zi)}, compute 



X 



J2 m-^i) Cov'i [F(z,), F{z,)\ AF{z,) (17) 



AF{z,) = F{z,)~FP{z,\s), 



(18) 



with FP{zi\s) = (dL(z|s)/Mpc)"^ 

For the sample of SNe we use in this study, 
we flux-averaged the SNe with dz — 0.07, to 
ensure that all redshift bins contain at least one 
SN. Our SN flux-averaging code is available at 
[http:/ /www.nhn.ou.edu/^wang/SNcode/, 



B. Galaxy Clustering Data 

For GC data, we use the distance measurements from 
the SDSS DR7 data. We compare two sets of distance 
measurements from SDSS DR7: (1) "CW2": the mea- 
surement oi H{z)rs{zd) and rs{zd)/ Da{z) (where H{z) is 
the Hubble parameter, Da{z) is the angular diameter dis- 
tance, and rs{zd) is the sound horizon at the drag epoch) 
by Chuang & Wang (2011) '28'] from the two-dimensional 
two-point correlation function of SDSS DR7 Luminous 
Red Galaxies (LRGs); (2) "WP2": the measurement of 
rs{zd)/Dv{z = 0.2) and rg{zd)/ Dv{z = 0.35) (where 
Dv{z) is the spherically-averaged distance) by Percival 
et al. (2010) [23| from the spherically-averaged power 
spectrum of combined data of SDSS DR7 LRG and main 
galaxy samples and 2-degree Field Galaxy Redshift Sur- 
vey (2dFGRS). 

Using the two-dimensional two-point correlation func- 
tion of SDSS DR7 in the scale range of 40-120 Mpc//i, 
Chuang & Wang (2011) ^ found that 

H{z ^Q.ib)rs{zd) = 13020±530km/s (4%) 
rs{zd)/DA{z = 0.i^) = 0.1518 ±0.0062 (4%) 

r = -0.0584, (19) 

where r is the normalized correlation coefficient be- 
tween H{z = 0.35)rs(zd) and rs{zd)lDA{z = 0.35), and 
rs{zd) is the sound horizon at the drag epoch (given 
by Eqs.ipS]) and ([5^ . The inverse covariance matrix of 
H{z — 0.35)rs(Z(i) and rs{zd)/ Da{z = 0.35) is 



GC,CW2 



0.35722E:- 
0.17833E: 



05 0.17833E: 
01 0.26104E:- 



01 
05 



(20) 



Spherically-averaged data give a measurement of dz 

rs{zd)/Dv{z), with 



Dviz) 



r{z)'^ cz 



1/3 



(21) 



Hiz) 

Chuang & Wang (2011) found that 

d^Z = rs{zd)/Dv{z = 0.35) = 0.1161 ± 0.0034 (2.9%) 

. (22) 

from the spherically-averaged correlation function of 
SDSS LRGs. 
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Using the spherically-averaged power spectrum of com- 
bined data of SDSS DR7 LRG and main galaxy samples 
and 2dFGRS, Percival et al. (2010) ^ found that 



,WP 
"0.2 

"0.35 



r,{zd)/Dviz = 0.2) = 0.1905 ± 0.0061 (3.2%) 
r^{zd)/Dv{z = 0.35) = 0.1097 ± 0.0036 (3.3%) 
0.337, (23) 



where r is the normalized correlation coefficient of do. 2 
and (io.35- The inverse covariance matrix of (do. 2, do. 35) 

is m 



"-"GCWP 



30124 -17227 
-17227 86977 



(24) 



For comparison, we also consider the distance measure- 
ments from the WiggleZ survey at z = 0.6 by Blake et 
al. (2011) and from the 6dF GRS at z = 0.106 by 
Beutler et al. (2011) 

das = r,(zd)/7:'y(z = 0.6) = 0.0692 ±0.0033 
do.106 = r,(zd)/i:'y(z = 0.106) = 0.336 ±0.015(25) 

GC data are included in our analysis by adding the 
following term to the of a- given model with 

xhc = [CgUp^^Pj)] Ap„ Ak = K-pf (26) 

For the Chuang & Wang (2011) GC measurements [28j . 
Pi = H{z = 0.35)rs{zd) and p2 = rs{za)/DA{z = 0.35). 
For the Percival et al. (2010) GC measurements |29| , 
Pi = do. 2 and p2 = do. 35. Note that pf°-*'°- (i = 1,2) are 
given in Egs. p^ and (|23p . and the inverse covariance 
matrices are given in Egs. ipO)) and ([M)) . Since the Wig- 
gleZ and 6dF surveys give independent measurements, 
Eg.d^ is replaced with Xgc = Yl^i=ii^Pi)'^ / ^pi^ ^\t\\ 



Pi = do. 6, and p2 
Eq.dlSI). 



^0.106: 



and CTp; (i = 1,2) 



Pi 

given in 



C. CMB data 

CMB data give us the comoving distance to the 
photon-decoupling surface r(z*), and the comoving sound 
horizon at photo-decoupling epoch rs{z^) [12] ■ Wang & 
Mukherjee 2007 [Hi showed that the CMB shift param- 
eters 



R = ^nrnHir{z,)/c, 
la = 7rr(z,)/rs(z*), 



(27) 



together with LUf, = ftbh^, provide an efficient summary of 
CMB data as far as dark energy constraints go. This has 
been verified by [3^ . Replacing Wf, with z» gives identical 
constraints when the CMB distance priors are combined 
with other data 1351. 



The comoving sound horizon at redshift z is given by 







a 



E{z'Y 



da' 



(28) 



+ Rba')a'''E-^{z') 



where a is the cosmic scale factor, a = 1/(1 + z), 
and a^i?^(z) = ri„i(a + Ocq) + ilfeC^ + VLxX{z)a'^, 

with Acq = ilrad/f^m = 1/(1 + ^oq), and Zoq = 

2.5 X lQ'^nmh'^{TcMB/'2-T^y^- The sound speed is 



Cs = 1/ y 3(1 + Rb a), with i?f, a = 3pf,/(4p^), Rb = 
SlbOOnbh^iTcMB/'Z.TK)-'^. We take Tcmb = 2.725. 

The redshift to the photon-decoupling surface, z,, is 
given by the fitting formula [11] : 



1048 [1 + 0.00124(f]b/i2)-o.738j _^ g^{n^hy^] , 

(29) 



where 



91 

92 



0.0783 {ftbh 



2\-0.238 



1 + 39.5 (f^6/l2)0.763 

0.560 

l + 21.1(f^6/l2)l-81 



(30) 
(31) 



The redshift of the drag epoch Zd is well approximated 
by M 



Zd 



where 



1291(r2™/i2)C 



1 + 0.659(f7™/i2)o.i 



[l + b^{^lbh^f^], (32) 



hi = 0.313(a„/^')""-^'Ml + 0.607(an/i')°-^^1 ,(33) 
&2 - 0.238{n,nhy-^^^. (34) 

We have derived the one-dimensional marginalized 
probability distributions (pdf) of (la, R, z^:,ns,rs{zd)) 
from WMAP7 data, for three different assumptions 
about dark energy: (1) wx{z) = —1; (2) wxiz) — w 
(constant); (3) wxiz) — wq + Wa{l — a) [1^. We find 
that these pdf 's are nearly independent of the assump- 
tion about dark energy, and are well fitted by Gaussian 
distributions with the following means and standard de- 
viations: 

(la) = 302.35, cr(L) = 0.86 
(R) = 1.728, cr(i?) = 0.02 
(z*) = 1091.32, ct(zO = 1.0 
{us) = 0.963, cr(n,) = 0.016 

{rs{zd)) = 152.85 Mpc,a(r,[zd]) - 1.85 Mpc. (35) 

Our Gaussian fits are fully consistent with the bestfit 
values for a constant w from Komatsu's website 1391. 
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The normalized 



covariance 



matrix 



of 



/ 1.0000 .20794 .47422 -.54889 .34914 \ 
.20794 1.0000 .74409 -.44370 -.76929 * 
.47422 .74409 1.0000 -.79575 -.19121 
-.54889 -.44370 -.79575 1.0000 -.09973 
V .34914 -.76929 -.19121 -.09973 1.0000 / 

(36) 

We find tliat the addition of rs{zd) is redundant, and 
amounts to partially double-weighting the CMB con- 
straints; so its measurement cannot be used in the CMB 
priors unless it is used to replace In(z^). This was not ap- 
parent when WMAP5 data were used Hoj; this is because 
WMAP7 data provide significantly tighter constraints on 
the set of CMB distance priors. 

Since the primary GC data we use in this paper have 
been marginalized over Us [28l |. we should marginalized 
the CMB distance priors over Ug as well. When marginal- 
ized over Us and rs{zd)^ the inverse covariance matrix of 
ila,R,z^) from WMAP7 is 



Coy 



-1 

CMB 



1.8571 
25.929 
-1.1433 



25.929 
5963.3 

-99.319 



-1.1433 
-99.319 
2.9443 



(37) 



CMB data are included in our analysis by adding the 
following term to the of ^ given model with pi = 
Liz^,), P2 = i?(z*),and = z*: 



XCMB 



Apj [CoV(jIj g {p., ,p-i)] Apj , Ap, = -p[ 



data 
i ' 

(38) 

where pf'^^'^ are the mean from Eq. (|35l) , and Cov^^^^ is 
the inverse covariance matrix of [i?(z,), Za(2*), z,] from 
Eq. ([37| . Note that p4 — Ug should be added if the con- 
straints on Us are included in the GC data. Finally, our 
Gaussian fit for the pdf of rs{zd) should be useful when 
fitting for BAO peak locations. 



D. Gammay-ray Burst Data 

We add gammay-ray burst (GRB) data to our analy- 
sis, since these are complementary in redshift range to the 
SN la data. We use GRB data in the form of the model- 
independent GRB distance measurements from Wang 
(2008c) [111, which were derived from the data of 69 
GRBs with 0.17 < z < 6.6 from Schaefer (2007) 

The GRB distance measurements are given in terms of 

sa 



oiZt) 



rp(0.17)' 



rp{z) 



(1 + z)i/2 i7o 



ch 



r(z). 



where r(z) is the comoving distance at z. 

The GRB data are included in our analysis by adding 
the following term to the of a given model: 



xIrb = [Arp(zi)] • (Covq^^) . • [Arpizj)] 



Ar^{z,) = 7v^^'^(z,) - r>(zO 



(40) 



where rp(z) is defined by Eq. (j39p . The covariance matrix 
is given by 

{CovGRB)ij = cr{rY,{z,))air^{zj)) (Cotgab)^^. , (41) 

where Coygrb is the normalized covariance matrix from 
Table 3 of Wang (2008c) [HI, and 

a(rv(zO) = <y (r>(z,0)+ , if rv(z) > r^izf^'^; 
fT(7v(zO) = a(j^{z,)r, if 7v(^) <7v(z)'^-'^(42) 

where a (f^(zi))^ and a (r^{zi))~ are the 68% C.L. errors 
given in Table 2 of Wang (2008c) HH. 



E. Dark energy parametrization 

Since we are ignorant of the true nature of dark energy, 
it is useful to measure the dark energy density function 
X{z) = px{z)/px{0) as a free function of redshift [4|3- 
|49| . This has the advantage of allowing dark energy mod- 
els in which px{z) becomes negative in the future, e.g., 
the "Big Crunch" models [50, 51], which are precluded 
if we parametrize dark energy with an equation of state 
wx{z) M- 

Here we parametrize X{z) by cubic-splining its values 
at z = 1/3, 2/3, and 1.0, and assume that X{z > 1) = 
X(z = 1). For simplicity of notation, we define X0.33 = 
X{z = 1/3), Xo.67 = X{z = 2/3), and Xi,o = X{z = 
1). Fixing X{z > 1) refiects the limit of current data, 
and avoids making assumptions about early dark energy 
that can be propagated into artificial constraints on dark 
energy at low z (sj, HI] . 

For comparison with the work of others, we also con- 
sider a dark energy equation of state linear in the cosmic 
scale factor a, wx(a) = wq + {1 — a)wa [11]. A related 
parametrization is [35| 



wx{z) = wo{3a - 2) + 3wq,5{1 - a). 



(43) 



(39) 



where wo.5 = wx{z = 0.5). Wang (2008b) [33 showed 
that (wo, Wo. 5) are much less correlated than (wq, Wq), 
thus are a better set of parameters to use. We find 
that (wo, Wo. 5) converge much faster than (wo, Wa) in 
a Markov Chain Monte Carlo (MCMC) likelihood anal- 
ysis for the same data. 



III. RESULTS 



^ The proper calibration of GRBs is an active area of research. 
For recent studies on the impact of detector thresholds, spectral 



analysis, and unknown selection effects, see, e.g., |43h4i 



We perform a MCMC likelihood analysis [52[ to obtain 
O(IO^) samples for each set of results presented in this 
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■ flux-averaged; stat+sys 

not flux-averaged; stat+sys 
- not fiux-averaged; stat only 




■ flux-averaged; stat+sys 
not flux-averaged; stat+sys 
not flux-averaged; stat only 

1 23.8 23.9 24 
M 



CMB+H +GRB+GC(CW2)+ 




— SNe flux-averaged 
SNe not flux-averaged 



CMB+H +GRB+GC{CW2)+ 




—SNe flux-averaged 
SNe not flux-averaged 



FIG. 1: The 2D marginalized contours of {w, Qm, yV() for SNe data 
(with and without fiux-averaging) , assuming a flat universe. The 
contours are at 68% and 95% confidence levels. 



FIG. 2: The 2D marginalized contours of (wa,Wa) and (uiqjW'o.s) 
for SNe data (with and without flux-averaging) combined with 
galaxy clustering (CW2), CMB, i^o, and GRB data. The contours 
are at 68% and 95% confldcnce levels. 



paper. We assume flat priors for all the parameters, and 
allow ranges of the parameters wide enough such that fur- 
ther increasing the allowed ranges has no impact on the 
results. The chains typically have worst e- values, defined 
to be the variance(mean)/mean(variance) of 1/2 chains, 
much smaller than 0.005, indicating convergence. The 
chains are subsequently appropriately thinned to ensure 
independent samples. 

In addition to the SN la, CMB, GC, and GRB data 
discussed in SecHIl we impose a prior of Hq = 73.8 ± 
2.4kms~^Mpc~"'^, from the HST measurements by Riess 
et al. (2011) [H. 

We do not assume a flat universe unless specifically 
noted. In addition to the dark energy parameters de- 
scribed in Sec lIIEl we also constrain cosmological pa- 
rameters (r^m, ^^fc, ^, Wb), where Wfc = ilfc/i^. In addi- 
tion, we marginalize over the SN la nuisance parameters 

We will present results for dark energy density at z — 
1/3, 2/3, and 1, as well as {wo,Wa) and {wo,wo,5), and a 
constant dark energy equation of state w. 



CMB+H +GRB+GC(CW2)+ 




jt riux-avaraged 

0.6 0.8 1 1.2 



CMB+H +GRB+GG|CW2)+ 
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-4 
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FIG. 3: The 2D marginalized contours of 

(+f 0.33, +f 0.67) +'^1.0, ^^m, f!fc) for SNe data (with and without 
flux-averaging) combined with galaxy clustering (CW2), CMB, 
Hq, and GRB data. The contours are at 68% and 95% confldence 
levels. 



A. The effect of flux-averaging of SNe la 

FigsH] shows the 2D marginalized contours of 
(w, rim,A4), assuming a constant equation of state for 
dark energy, w, and a flat universe. Note that the 
inclusion of systematic errors of SNe leads to signifi- 
cantly larger uncertainties in estimated parameters, com- 
pared to when only statistical errors of SNe are included. 
Clearly, flux-averaging (thick solid lines) leads to larger 
errors on dark energy and cosmological parameters if only 
SN la data are used. 

Fig[2] shows the 2D marginalized contours of {wo,Wa) 
and {wo,wo.5) for SNe data combined with with CMB, 
Ho, GRB, and GC (CW2) data. The solid contours 
are for flux-averaged SNe, and dotted contours are for 
SNe without flux-averaging. The SNe data with or with- 
out flux-averaging lead to qualitatively consistent results, 
with flux-averaging expanding the parameter space at 
wq < —1, which shifts the bestfit model from wq > —1 
towards wq = —1, and from Wa < towards Wa — (it 



has a smaller impact on wq.s since wq.s is less correlated 
with Wo than Wa). Table U tabulates the mean and 68.3% 
confidence intervals for (wojWa) and {wojWo,^) for SNe 
data combined with with CMB, Hq, GRB, with or with- 
out GC (CW2) data (the latter corresponds to Fig[21). 

FigO shows the 2D marginalized contours of 
(Xo.33, +Yo.67, ^i.Oj f^m, f^fc) for SNc data combined with 
CMB, Ho, GRB, and GC data (CW2). The solid con- 
tours are for fiux-averaged SNe, and dotted contours 
are for SNe without fiux-averaging. The SNe data with 
or without flux-averaging lead to qualitatively consis- 
tent results, with a shift of the bestfit model towards 
a cosmological constant; this is the same trend as in the 
{wq, Wa) and (wq, wq.s) parametrizations. Note that flux- 
averaging leads to signiflcantly tighter constraints on the 
dark energy density at z = 1, Xi.o; this indicates that 
the impact of flux-averaging increases with redshift, since 
(+^^^0.337 +^o.67j -''^i.o) are only weakly correlated. This is as 
expected, since the bias in estimated parameters due to 
weak lensing increases with redshift 27 1. 
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TABLE I: Effect of flux-averaging SNe on (wo, Wa) and (wojWo.s) 



flux-averaging wo Wa wo wo.s 

SNe-HCMB-l-iyo+GRB 

yes -0.987 (-1.247, -0.727) -0.780 (-2.046, 0.516) -0.997 (-1.268, -0.722) -1.260 (-1.546, -0.965) 

no -0.780 (-1.013, -0.545) -1.424 (-2.891, 0.045) -0.783 (-1.023, -0.547) -1.248 (-1.549, -0.946) 

SNe-hCMB+i7o+GRB-hGC(CW2) 

yes -0.995 (-1.211,-0.776) -0.676 (-1.884, 0.566) -1.000 (-1.216, -0.784) -1.211 (-1.454, -0.966) 

no -0.832 (-1.061,-0.606) -1.353 (-2.844, 0.151) -0.842 (-1.067, -0.615) -1.266 (-1.559, -0.972) 



B. Comparison of different galaxy clustering data 

We now explore the differences of the various galaxy 
clustering (GC) measurements in constraining dark en- 
ergy, assuming a constant equation of state for dark 
energy, w. FigU] shows the 2D marginalized contours 
of (w,r2m,f2fc) for different GC measurements combined 
with CMB, Ho, and GRB data. 

The first row of FigU compares the H(z ~ 0.35)rs{zd) 
and Tsizd)/ Da(z = 0.35) measurements by Chuang & 
Wang (2011) Q with their do.35 = rs{zd)/Dv{z = 0.35) 
measurement (both from SDSS DR7 LRGs), as well as 
the do. 2 and do. 35 measurements by Percival et al. (2010) 
[29} from SDSS DR7 LRG and main galaxy samples 
and 2dFGRS, and the do. 6 measurement by Blake et al. 
(2011) from the WiggleZ survey [13] combined with the 
do. 106 measurement by Beutler et al. (2011) from 6dF 
GRS. 

For the Chuang & Wang (2011) ^ GC measure- 
ments (CW2 and CWl), the constraints on w are tight- 
ened significantly by going from spherically-averaged 
data (CWl), i.e., do.35, to 2D data (CW2), i.e., H{z = 
0.35)rs{zd) and rs{zd)/ Da{z = 0.35), as indicated by 
comparing the thin solid contours (CWl) to thick solid 
contours (CW2) in the first row of Fig|4l This is as ex- 
pected, as more information from GC is included in CW2 
compared to CWl. 

Both the Percival et al. (2010) GC measurements 
(WP) and the combined WiggleZ survey and 6dF GRS 
measurements (CB-|-) favor w < —1, while the Chuang 
& Wang (2011) [11] GC measurements favor w — —1. 

The second row in Fig|4] compares the do. 2 and do.35 
measurements by Percival et al. (2010) ^ (WP2), with 
their measurements of do. 2 and do.35 separately. Clearly, 
most of the constraining power on w comes from do.35. 
While the do. 2 measurement favors w = —I, the do.35 
measurement favors w < —I. 

The measurements of do.35 by Chuang & Wang (2011) 
[2^ and Percival et al. (2010) ^29] are similar in preci- 
sion, but differ systematically: d^^ = rs{zd)/ Dv{z ~ 
0.35) = 0.1161 ± 0.0034, while df£ = r,{zd)/Dvlz = 
0.35) = 0.1097 ± 0.0036. The lower measured value of 
implies a smaller H{z — 0.35), which in turn im- 
plies a more negative w. When combined with CMB, Hq, 
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FIG. 4: The 2D marginalized contours of (ui, Hj.) for differ- 
ent galaxy clustering measurements combined with CMB, Hq, and 
GRB data. The contours are at 68% and 95% confidence levels. 



and GRB data, favors w 



135 idvuio uj — —1, while d^3^ favors 
w < ~1. Note that these two measurements used differ- 
ent methods to analyze GC data: Chuang & Wang (2011) 
used the galaxy correlation function, while Percival et al. 
(2010) used galaxy power spectrum. It is not surprising 
that they lead to different distance measurements from 
GC. 

For the rest of this paper, we will present only results 
using the GC measurements by Chuang & Wang (2011) 
[lH; these are more conservative as they were derived 
without assuming CMB priors. 



C. Constraints on dark energy and H{z) 

Figs l5][7] show the 2D marginalized contours of dark en- 
ergy and cosmological parameters for the three different 
sets of dark energy parameters: (1) wx{z) = u'=const.; 
(2) wxiz) = Wo + Wa(l — a), and wxiz) — ?x;o(3a — 2) + 
3^0.5(1 - a); (3) X{z) = px{z)/px{0) given by a cu- 
bic spline at z < 1 over X^. = X{zi), with zi — 1/3, 
Z2 = 2/3, and Zg = 1, and X(z > 1) = X{z = 1). 

FiglS] shows the dark energy density function px(z) 
measured from SNe data (flux-averaged with dz = 0.07), 
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FIG. 5: The 2D marginalized contours of (ui, flm, ^k) for SNe data 
(flux-averaged), galaxy clustering measurements (CW2), combined 
with CMB, Ho, and GRB data. The contours are at 68% and 95% 
confldcnce levels. 



CMB+H +GRB+ 




-SNe(t.a.)+GC(CW2) 
-SNe(t.a.) 
GG(CW2) 



CMB+H^+GRB+ 




-SNe(f.a.)+GG(CW2) 
-SNe(f.a.) 
GC(CW2) 



SNe+GC+CMB+Ho+GRB 

Itijao 





0.2 0.4 o.a 0.8 



FIG. 6: The 2D marginalized contours of {wo,Wa) and (uiOjtoo.s) 
for SNe data (flux-averaged), galaxy clustering measurements 
(CW2), combined with CMB, Ho, and GRB data. The contours 
are at 68% and 95% confidence levels. 



galaxy clustering measurements (CW2), combined with 
CMB, Ho, and GRB data. The upper panel of FigE] 
shows both px{z) and pm{z), in units of yoctograms per 
cubic meter, with 1 yoctogram = 10~^^ grams. The un- 
certainties in X{z), VLx = 1 — — ^k, and h have been 
propagated into that of px{z)- The lower panel of FigE] 
shows the corresponding X{z) — px{z)/px{0)- FigE] 
shows the corresponding cosmic expansion history H{z). 
Clearly, given the same priors of CMB, Hq, and GRB 



Px(z) 
Px(0) 





0.1 0.2 0.2 0.25 0.3 0.2 0.25 0.3 

n, n n 

k m m 

FIG. 7: The 2D marginalized contours of 

(^0.33i ^0.67i ^l.Oi f^m, f!fc) for SNe data (flux-averaged), 
galaxy clustering measurements (CW2), combined with CMB, 
Ho, and GRB data. The contours are at 68% and 95% confidence 
levels. 



FIG. 8: Dark energy density function px{^) measured from SNe 
data (fiux-averaged with dz = 0.07), galaxy clustering measure- 
ments (CW2), combined with CMB, Ho, and GRB data. The 
contours are at 68% and 95% confidence levels. The upper panel 
shows both px{^) S'lid pm{z) in units of yoctograms per cubic me- 
ter, with 1 yoctogram = 10"^* grams. The lower panel shows the 
corresponding X{z) = px{z)/ Px{0)- 



data, SNe lead to much stronger constraints on dark en- 
ergy than galaxy clustering data at present. Note also 
that the addition of galaxy clustering data to SN data 
leads to significantly improved constraints. 

Table HI] gives the dark energy density function X{z) = 
px{z)/ px{0) and the cosmic expansion history H{z), 
as well as {Q.m,^k,h,u}ii), measured from current data 
of SNe-t-GC-hCMB+Ho+GRB. The H[z) measurements 
are derived using Eq.([21). Tables illllllVI give the normal- 
ized covariance matrices of the X{z) and H{z) measure- 
ments. Note that both the X{z) and H{z) measurements 
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FIG. 9: The cosmic expansion history H(z) measured from SNe 
data (flux-averaged with dz = 0.07), galaxy clustering measure- 
ments (CW2), combined with CMB, Hq, and GRB data. The 
error bars indicate the 68.3% confidence level ranges. 



are only weakly correlated at different redshifts. 

To quantify the comparison in dark energy constraints, 
we can use the general dark energy Figure-of-Merit 
(FoM) defined by Wang (2008b) [H 



FoM = 



1 



v/detCov(/i,/2,...,/Ar) 



(44) 



where (/i, /2, /jv) is the set of parameters that have 
been chosen to parametrize dark energy. The Dark En- 
ergy Task Force (DETF) defined the dark energy FoM 
to be the inverse of the area enclosed by the 95% confi- 
dence level contour of {wo,Wa) [I^. The areas enclosed 
by contours are difficult to calculate for real data, as 
these contours can be quite irregular. The definition of 
Eq. (|44l) has the advantage of being easy to calculate for 
either real or simulated data. For (/i,/2) = {wo,Wa), 
FoM of Ea. (l44l) is proportional to the FoM defined by 
the DETF for ideal Gaussian-distributed data, and the 
same as the relative FoM used by the DETF in Fisher 
matrix forecasts. 

Table |V] shows the dark energy FoM from SNe (fiux- 
averaged), galaxy clustering measurements (CW2), to- 
gether with CMB, Hq, and GRB data, for (wo,Wa), 
(wcwo.s), and (X0.33, Xo.67, -'^i.o)- 



IV. SUMMARY AND DISCUSSION 

We have examined how dark energy constraints from 
current observational data depend on the analysis meth- 
ods used: the analysis of Type la supernovae (SNe la), 
and that of galaxy clustering data (GC). 



We generalize the flux-averaging analysis method of 
SNe la to allow correlated errors of SNe la, in order 
to reduce the systematic bias due to weak lensing of 
SNe la. We find that fiux-averaging leads to larger 
errors on dark energy and cosmological parameters if 
only SN la data are used (see Fig[T]). When SN la 
data (the latest compilation by the SNLS team) are 
combined with WMAP 7 year results, the latest Hub- 
ble constant (Hq) measurement using the HST, and 
gamma ray burst (GRB) data, fiux-averaging of SNe la 
increases the concordance of SNe la with other data, and 
shifts the bestfit cosmological model notably closer to 
w — —1 (see Fig [2] and Table This leads to signifi- 
cantly tighter constraints on the dark energy density at 
z = 1, and the cosmic curvature ftk (see FiglS]). We have 
made our SN flux-averaging code publicly available at 



http: / / www.nhn.ou.edu/ ~wang/SNcode/ 



Note that since the flux-averaging of SNe la increases 
the concordance of SN la with other data, the combined 
data with flux-averaging of SNe la gives comparable con- 
straints for most parameters, and smaller uncertainties 
for dark energy density at z = 1, Xi o (most affected by 
weak lensing), and fife (strongly correlated with Xi.o), 
compared with combined data with no flux-averaging of 
SNe la. The combination of concordant data tightens 
the overall constraints, while that of discordant data may 
not. 

We find that given the same pirors (CMB+iJo+GRB) 
both the Percival et al. (2010) GC measurements |29j 
(WP: do. 2 and do. 35) and the combined WiggleZ survey 
and 6dF GRS measurements [13, |3l| (CB-|-: do.e and 
1^0. loe) favor w < —1^, while the Chuang & Wang (2011) 
GC measurements [28| (CW2: H{z = 0.35)r,(zd) and 
rs{zd)/DA{z = 0.35)) favor w = -1 (see Figg]). This 
could be due to the difference in the analysis methods 
used to derive these GC measurements. We find that the 
GC measurements by Chuang & Wang (2011) are consis- 
tent with SN la data (which favor w — —1), and tighten 
the dark energy constraints significantly when combined 
with SN data (see FigEHTl). 

For the convenience of future work, we have derived 
Gaussian fits to the probability distributions of a set of 
CMB parameters, [la{z*), R{z^,), z^,ns,rs{zd)], that sum- 
marize the WMAP7 data, as well as well their covari- 
ance matrix. Note that while [la{z*), R{z■^,), z^] are suf- 
ficient for summarizing the CMB priors when Ug de- 
pendence has been marginalized over in GC analysis, 
[laizt:), R{z^:), Zf,,ns] need to be used when the con- 
straints on ris from GC are included (as will be the case 
for future GC data). We find that rs{zd) gives redun- 
dant information, its measurement cannot be used in the 
CMB priors unless it is used to replace la{z*)- Also, our 
Gaussian fit for the pdf of rs{zd) from WMAP7 can be 



^ Similar results were found by [54 
GC measurements [2911. 



[sa l using Percival et al. (2010) 
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TABLE II: X{z), H{z), and cosmological parameters estimated from current data with 68.3% C.L. upper and lower uncertain- 
ties. 







a 


fi + a 


ji — (j+ 


X{z = 


1/3) 


0.969 


0.108 


0.862 


1.077 


X{z = 


2/3) 


1.152 


0.347 


0.812 


1.492 


X{z = 


1.0) 


0.453 


0.863 


-0.385 


1.284 






0.252 


0.016 


0.236 


0.268 






-0.0049 


0.0131 


-0.0180 


0.0080 


h 




0.734 


0.019 


0.715 


0.753 


UJb 




0.02234 


0.00062 


0.02170 


0.02290 


H{z = 


l/i)/Ho 


1.148 


0.040 


1.107 


1.186 


H{z = 


2/S)/Ho 


1.419 


0.093 


1.328 


1.510 


H{z = 


1.0) /Ho 


1.511 


0.221 


1.289 


1.731 



TABLE III: Normalized covariance matrix of X{z) from current data 





1 


2 


3 


1 


1.0000 


0.1830 


-0.1195 


2 


0.1830 


1.0000 


-0.1283 


3 


-0.1195 


-0.1283 


1.0000 



useful when fitting for BAO peaks in GC data analysis. 

We find that current data are fully consistent with a 
cosmological constant and a flat universe, consistent with 
the latest findings by others (see, e.g., [l8l.[55l-t61j). Since 
the uncertainties remain large, a deviation from a cosmo- 
logical constant is far from being ruled out (see FigH]). 

At present, SN data provide much stronger constraints 
than GC data (see Table W} . This will change when am- 
bitious galaxy surveys are carried out from space in the 
future [62, 63]. We can expect dramatic progress in our 
measurement and understanding of dark energy in the 
next decade or so. 
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